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摘要 : FSET. RP PA AR OG EEE, 
视 为 未 来 核 聚 变 装置 中 最 有 前 景 的 面向 等 离子 材料 。 在 聚变 服役 环境 下 ，14Mey 的 高 能 中 
子 以 及 低能 氨 / 氨 粒子 流 对 钨 基 材 料 造 成 严重 的 辐 照 损伤 。 研 究 材料 的 辐 照 损伤 与 氧 氛 效 应 
机 理 对 揭示 辐 照 引起 材料 微观 结构 与 性 能 的 变化 以 及 探索 开发 新 型 抗 辐 照 材 料 具有 重要 的 
意义 。 近 年 来 ， 随 着 计算 模拟 技术 的 发 展 ， 多 尺度 模拟 方法 在 聚变 堆 关 键 材料 辐 照 损伤 与 
氧 氨 效 应 机 理 研究 方面 有 着 广泛 的 应 用 。 本 文 主要 结合 作者 近 几 年 的 研究 实践 ， 介 绍 了 纳 
米 晶 钨 / 铁 的 抗 辐 照 损 伤 机 理 、 钨 中 氢 滞 留 机 制 、 合 金 元 素 强化 界面 机 制 、 妨 变 元 素 偏 聚 及 
其 对 氢 泡 形成 机 制 的 影响 。 通 过 多 尺度 模拟 方法 来 获取 材料 在 辐 照 环境 下 的 不 同时 间 / 空 间 
层次 的 结构 演化 信息 ， 从 而 揭示 材料 的 辐 照 损伤 与 氨氮 效应 机 理 ， 为 聚变 堆 面向 等 离子 体 
材料 的 优化 设计 与 性 能 调控 提供 理论 基础 。 
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一 、 引 言 

随 着 传统 化 石 能 源 的 逐渐 枯竭 以 及 人 类 对 能 源 需求 的 不 断 增长 ， 能 源 问 题 已 成 为 人 类 
生存 与 发 展 的 重大 问题 之 一 。 核 聚变 能 被 认为 是 可 以 最 终 解决 人 类 能 源 问题 的 重要 途径 之 
一 。 利 用 强 磁 场 约束 高 温 等 离子 体 的 托 克 马克 《Tokamak) 是 最 有 希望 实现 受 控 热 核 聚变 
反应 的 装置 [1]。 由 于 其 技术 复杂 和 条 件 苛刻 ， 磁 约束 托 卡 马克 装置 仍 面 临 一 些 急 需 解决 的 
(S) 关键 问题 。 其 中 ， 聚 变 堆 材料 问题 ， 尤 其 是 面向 等 离子 体 材料 (PFMs)， 是 制约 托 卡 马克 装 
置 性 能 及 其 发 展 的 关键 问题 之 一 [2]。PFMs 作为 直接 面 对 高 温 等 离子 体 的 第 一 壁 、 偏 滤器 
等 的 护 甲 材 料 ， 其 工作 环境 极其 苛刻 ， 遭 受 着 高 温 、 高 热 负 荷 、 强 束 流 粒子 与 中 子 辐 照 等 
综合 作用 。 研 究 表明 国际 热 核 聚变 实验 堆 CTER) 偏 滤 器 要 承受 极 高 的 热流 〈10-20 
MW/m”)， 以 及 强 束 流 (10”-10* m?s) 和 低能 (<100 eV) 离子 流 的 辐 照 [3]。 

E CW) 以 其 高 熔点 、 低 溅 射 率 和 高 热 导 率 等 优点 而 被 视 为 未 来 聚变 堆 中 最 有 前 景 的 
PFM[4-5]。ITER 和 中 国 东方 超 环 (EAST， 中 科 院 等 离子 体 物 理 研究 所 ) 装置 中 已 经 使 用 
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纯 钨 作为 偏 滤器 材料 。 然 而 ， 在 聚变 服役 过 程 中 ， 高 能 中 子 将 对 钨 基 材 料 造 成 严重 辐 照 损 
伤 ， 导 致 材料 中 产生 大 量 的 辐 照 缺陷 〈 如 空位 和 自 间 隙 原子 ) 。 这 些 缺 陷 可 以 迁移 、 聚 集 
形成 间隙 团 、 堆 埃 四 面体 以 及 空洞 等 ， 进 而 引起 材料 宏观 性 能 改变 ， 如 肿胀 、 硬 化 、 脆 化 
以 及 非 唱 化 等 ， 最 终 导致 材料 服役 性 能 降低 甚至 失效 [6-9]。 此 外 ， 大 量 的 低能 高 通 量 氧 、 
所 粒子 的 辐 照 以 及 由 妨 变 反应 而 产生 的 大 量 氧 、 氢 也 会 对 材料 的 结构 与 性 能 造成 严重 影响 
在 聚变 服役 工 况 条 件 下 ， 钨 基 材 料 会 同时 受到 高 剂量 离 位 损伤 、 氢 氨 效 应 以 及 妨 变 效应 驴 
共同 作用 。 因 此 ， 材 料 的 辐 照 损 伤 效应 是 一 个 极其 复杂 的 过 程 。 材 料 的 微观 组 织 、 结 构 缺 
陷 、 辐 照 剂 量 、 氧 / 氧 聚 集 、 外 加 温度 等 多 种 因素 共同 决定 材料 的 结构 和 宏观 力学 性 能 。 单 
纯 依赖 现 有 的 实验 方法 和 检测 技术 ， 尚 难以 系统 而 全 面 地 认识 材料 微观 结构 与 宏观 性 能 之 
间 的 内 在 联系 ， 特 别 是 实验 上 难以 直接 跟踪 和 探测 原子 尺度 下 缺陷 的 演化 行为 。 借 助 现代 
先进 的 计算 模拟 技术 来 研究 材料 的 辐 照 损伤 微观 机 理 已 经 成 为 男 外 一 条 重要 的 途径 。 多 尺 


到 度 模拟 技术 已 广泛 应 用 于 钨 基 材料 的 相关 研究 中 [10-12]。 其 中 : 原子 尺度 的 第 一 性 原理 方 
他 法 研究 材料 中 点 缺陷 /杂质 原子 的 形成 能 、 扩 散 路 径 与 激活 能 、 复 合 团 簇 的 稳定 性 、 聚 集 和 
c0 解 离 行 为 、 缺 陷 与 界面 间 相互 作用 等 性 质 ， 微 观 尺度 的 分 子 动力 学 模拟 初级 离 位 损伤 的 产 
S 生 、 级 联 位 移 过 程 、 点 缺陷 的 迁移 和 团聚 、 缺 陷 团 簇 的 稳定 性 和 迁移 率 、 杂 质 /缺陷 与 位 错 
© / 晶 界 之 间 的 相互 作用 等 ， 微 观 和 介 观 尺度 的 蒙特 卡 洛 方法 和 速率 理论 方法 研究 缺陷 和 位 错 
z WPR L/S RIN TES. MARRE: PARERE AA 
= 材料 微 结构 与 宏观 力学 性 能 的 对 应 关系 ， 宏 观 尺 度 的 连续 介质 力学 和 有 限 元 等 方法 评估 到 
N 变 堆 各 功能 模块 在 实际 运行 环境 中 的 表现 。 


本 文 主要 介绍 近年 来 我 们 课题 组 利用 多 尺度 模拟 方法 研究 铝 中 辐 照 损伤 与 氢 氨 行为 机 
理 的 一 些 进展 ， 如 纳米 唱 铝 / 铁 的 抗 辐 照 损 伤 机 理 、 钨 中 氢 滞 留 机 制 、 合 金 元 素 强化 界面 机 
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二 、 纳 米 唱 钨 / 铁 抗 辐 照 损伤 机 理 跨 尺 度 模拟 研究 

近 十 年 来 ， 人 们 在 实践 中 逐渐 认识 到 纳米 结构 材料 通常 具有 较 好 的 抗 辐 照 损伤 性 能 
[13-16]， 这 与 缺陷 阱 (表面 、 唱 界 等 ) 对 辐 照 缺陷 的 捕获 、 进 而 促进 其 复合 有 关 。 为 研究 
纳米 结构 金属 材料 中 辐 照 缺 陷 与 界面 跨 尺度 相互 作用 ， 我 们 自主 开发 了 一 套 多 尺度 模拟 程 
Fe (AL) 。 模 块 包括 分 子 动力 学 MD) 程序 ， 温 度 加 速 动 力学 (TAD) 程序 ， 分 子 静 力 
学 (MS) 程序 ， 实 体 动力 学 蒙特 卡 洛 (OKMC) FER, SEIN THERE MBH A Se 
洛 程序 和 探索 体系 势能 面 的 “弹性 带 ” CNEB) 程序 , “两 体 ” (dimer) 程序 以 及 辐 照 损伤 
结构 可 视 化 与 分 析 等 辅助 程序 。 基 于 上 述 多 尺度 模拟 程序 ， 我 们 研究 了 钨 、 铁 晶 界 附近 辐 
照 产生 的 自 间隙 和 空位 的 扩散 、 偏 了 请、 复合 过 程 ， 揭 示 了 晶 界 促进 的 缺陷 偏 聚 和 复合 ，; 
而 提高 材料 抗 辐 照 损伤 的 机 理 [17-21]。 


图 1. 辐 照 缺陷 演化 跨 斥 度 模拟 程序 模块 。 


Figure 1. Multi-scale modeling program for irradiation defect evolution. 


使 用 分 子 动力 学 和 静态 计算 方法 ， 考 察 了 唱 界 在 修复 纳米 结构 钨 中 辐 照 损伤 方面 扮演 
的 能 量 学 和 动力 学 角色 。 研 究 发 现 ， 在 辐 照 过 程 中 ， 相 对 于 空位 而 言 ， 间 阶 原 子 被 唱 界 优 
先 吸 收 ， 从 而 在 级 联 碰撞 冷却 后 形成 唱 界 附近 富 含 空 位 而 品 界 内 部 间隙 高 度 局 域 的 缺陷 结 
构 。 唱 界 通过 降低 其 附近 空位 / 间 队 的 形成 能 和 扩散 能 垒 而 作为 缺陷 阱 ， 特 别 是 间隙 在 唱 界 
附近 不 需要 克服 能 爸 即 可 迁 入 唱 界 。 优 先进 入 唱 界 内 的 间隙 可 与 附近 的 空位 以 低能 又 、 多 
原子 协同 参与 的 方式 复合 。 同 时 我 们 发 现 ， 唱 界 增强 空位 扩散 、 复 合 的 区 域 较为 有 限 〈1- 
1.5 纳米 ) ， 使 得 该 区 域 体积 分 数 较 小 〈 几 个 百分点 ) 。 这 表明 实验 上 提高 钨 基 材 料 抗 辐 照 
性 能 的 两 种 途径 : 细 化 钨 唱 粒 和 加 入 纳米 小 颗粒 ， 后 者 可 能 更 为 有 效 [21]。 
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T 针对 候选 结构 材料 - 铁 素 体 / 马 氏 体 钢 的 辐 照 损伤 问题 ， 我 们 使 用 分 子 动力 学 、 静 态 计 
Ô 算 和 温度 加 速 的 动力 学 方法 ， 考 察 了 品 界 附近 小 型 空位 团 簇 的 能 量 学 和 动力 学 行为 。 研 究 
N RH, HER SEA ARI EMMER KA AMENDE RMR E AE 
= 数 1-9) 的 阱 。 空 位 团 簇 在 蝇 界 附近 具有 较 高 的 活性 。 其 经 过 块 体 扩散 、 局 部 吸收 和 局 部 沿 
o0 晶 扩散 过 程 而 被 唱 界 吸收 [20]。 为 了 深入 考察 晶 界 对 辐 照 缺陷 产生 和 演化 的 影响 ， 我 们 使 


~N 用 静态 计算 系统 考察 了 钨 、 钼 、 铁 等 金属 中 部 分 唱 界 附近 空位 / 间 阶 的 扩散 复合 问题 ， 获 得 
了 蝇 界 与 缺陷 作用 的 一 系列 参数 : 缺陷 形成 能 ， 偏 聚 能 ， 扩 散 -复合 能 垒 ， 以 及 相应 的 作用 
范围 ， 这 为 使 用 高 级 模拟 工具 葛 定 了 基础 [19]。 


综合 分 子 动力 学 、 静 态 计算 和 实体 动力 学 蒙特 卡 洛 方法 ， 考 察 了 铁 纳米 唱 中 辐 照 点 缺 

隐 与 晶 界 跨 尺度 基本 相互 作用 [17,18]。 重 点 关注 块 体 和 晶 界 内 部 的 缺陷 扩散 、 复 合 过 程 及 

= 其 之 间 的 耦合 作用 。 研 究 发 现 铁 纳米 唱 抗 辐 照 性 能 与 其 品 界 附近 和 内 部 耦合 的 辐 照 缺陷 仿 
© 聚 、 复 合 过 程 有 关 《〈 图 2) 。 当 自 间隙 原子 处 于 块 体 或 唱 界 附近 时 ， 在 其 周围 均 形成 空位 - 


自 间隙 低能 爸 复 合 区 域 ， 有 利于 空位 - 自 间隙 复合 。 复 合 涉及 两 个 基本 的 耦合 过 程 。 在 自 间 
隙 原子 和 空位 偏 聚 到 唱 界 前 ， 它 们 在 块 体 中 复合 了 相当 高 的 比例 。 未 在 块 体 中 复合 的 自 间 
隙 原子 偏 聚 至 晶 界 处 ， 以 沿 晶 扩 散 、 巡 游 的 方式 复合 晶 界 附近 的 空位 ， 而 不 必 发 射出 去 。 
这 两 种 机 理 均 能 在 低温 下 起 作用 。 该 工作 给 出 了 铁 纳米 晶 微观 和 宏观 的 辐 照 损伤 自修 复 图 
像 ， 如 图 3 所 示 。 


图 2. 分 子 动力 学 (molecular dynamics, MD) 模拟 1000K 时 铁 并 5(310X[001] 倾 侧 对 称 晶 界 
附近 辐 照 缺 陷 产 生 与 演化 过 程 。(a) 和 (b) 给 出 了 初级 碰撞 原子 PKA) 能 量 3keV， 分 别 距 
唱 界 32 和 23A 的 情形 。(c) 给 出 了 辐 照 产 生 的 间隙 (SIA〉、 空 位 (V) 与 晶 界 基 本 相互 作用 过 
程 。 包 括 间隙 空位 在 块 体 区 域 的 扩散 、 复 合 过程 1， 在 晶 界 附近 偏 聚 过 程 2 和 2’， 晶 界 附 

近 复 合 过程 3 和 3?， 以 及 在 晶 界 内 部 扩散 、 复 合 过 程 4。 刻 画 这些 过 程 的 参数 有 激活 能 (EJ 
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着 色 ， 颜 色 范围 从 蓝 到 红 对 应 势能 最 小 和 最 大 值 。 标 出 的 时 刻 时 间 单 位 为 皮 秒 。 


Figure 2. Representative simulation snapshots illustrating the progression of a collision cascade 
near a 5(310)[001] symmetric tilt Fe GB at 1000 K. In (a)-(b), the damage cascade was initiated 
by a single PKA with a kinetic energy of 3 keV placed at a distance of about 52 A and 23 A 
normal to the GB, respectively. (c) shows schematic illustration of the fundamental atomic 
interaction processes between the SIA, V and the GB, including diffusion and annihilation of the 
SIA and V in the bulk region far away from the GB (process 1), the segregation of the SIA and V 
near the GB (processes 2 and 2’, respectively), the SIA-V annihilation in the vicinity of the GB 
(processes 3 and 3’, respectively), and their migration and recombination along the GB (process 
4). These processes can be characterized according to their energetic and kinetic parameters such 
as activation energy barrier E., activation temperature 7,, segregation energy Eseg, and the range of 
interaction between the defects and the GB, W. In (a)-(b), the atoms are colored with their 
potential energy and the colors from blue to red correspond to the minimum of maximum of the 


potential energy. The time is given in ps. 


图 3. 结合 有 限 温 度 下 MD 模拟 ， OK 下 MS 计算 和 OKMC 模拟 计算 给 出 的 辐 照 缺陷 与 品 界 
作用 图 像 。 图 中 过 程 编 号 与 图 2(o) 一 致 。 绿 色 球 代表 间隙 ， 红 色 立 方 体 代表 空位 。 围 绕 绿 
色 球 的 红色 椭 球 代表 间隙 周围 形成 的 空位 -间隙 低能 垒 复合 区 域 。 


Figure 3. Schematic illustration of the interaction between irradiation defects and grain boundary 
in nano-crystalline iron by MD simulations at finite temperatures, MS calculations at 0 K and 
OKMC simulations at 0 K. Here, the numbers denote the same atomic processes as in Fig. 2c. The 
green sphere represents the SIA, the red cube denotes the V, and the orange sphere and ellipsoid 
around the green sphere represent the low-barrier annihilation region around an SIA in the bulk 


and at the GB, respectively. 
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低 氧 滞留 是 钨 作为 面向 等 离子 体 材料 的 必要 条 件 。 钨 中 氧 滞留 量 取决 于 材料 中 的 缺陷 
密度 ， 而 缺陷 密度 与 材料 加 工 和 服役 环境 密切 相关 [22]。 在 聚变 堆 环境 下 ， 辐 照 产生 的 大 
量 缺 陷 可 作为 氧 的 捕获 点 ， 致 使 钨 中 氧 滞留 量 升 高 ， 威 胁 聚 变 堆 的 安全 运行 。 实 验 研 究 表 
明 钨 辐 照 后 氧 滞留 行为 十 分 复杂 ， 与 材料 性 质 〈 杂 质 浓 度 、 合 金 元 素 种 类 、 唱 粒 尺 寸 、 热 
处 理工 艺 、 材 料 表面 性 能 等 ) 和 辐 照 条 件 〈 辐 照 预 处 理 过程 、 氧 同位 素 离子 辐 照 剂量 和 通 
EE 、 辐 照 温 度 等 ) 密切 相关 [23-29]。 但 是 由 于 实验 环境 和 条 件 的 差异 ， 实 验 结果 间 往 往 存 
在 很 大 的 差异 ， 甚 至 彼此 间 相 互 对 立 。 因 此 ， 和 急需 从 微观 角度 揭示 辐 照 环境 下 钨 基 材 料 所 
滞留 机 理 ， 理 解 实验 结果 间 的 差异 ， 进 而 全 面 认识 忽 基 材料 中 氧 滞留 行为 ， 为 降低 辐 照 环 


由 
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(DFT) 系统 研究 了 钨 中 间隙 氢 溶 解 、 扩 散 和 团 艇 化 行为 、 计 算 了 各 种 点 缺陷 〈 空 位 、 杂 
质 、 合 金 元 素 等 ) 与 所 相互 作用 参数 ， 建 立 了 完备 钨 中 氢 与 点 缺陷 相互 作用 数据 库 。 在 此 
基础 上 ， 我 们 进一步 开发 了 用 于 模拟 钨 中 氧 滞留 行为 的 OKMC 程序 ， 研 究 了 氧 离子 辐 照 后 
的 热 脱 附 行为 。 


图 4. 钨 中 氧 扩散 系数 模拟 和 实验 数据 比较 〈 图 中 Frauenfelder Zakharov, Benamati, Otsuka, 
Ikeda, Hoshihira, Heinola, Liu 和 Ikeda-Ostuka 为 实验 值 ，hTST+TE 和 WZ+TE 是 基于 过 渡 态 
理论 和 Wert-Zener 模型 并 考虑 温度 效应 而 得 到 的 理论 计算 值 ) [30]。 


Figure 4. Calculated diffusivity of hydrogen in tungsten in comparison with experimental data: (a) 
diffusion coefficient vs. temperature; (b) diffusion coefficient vs. the reciprocal of temperature; (c) 
diffusion coefficient vs. the reciprocal of temperature for the high temperature range. 
Experimental data are from Frauenfelder, Zakharov, Benamati, Otsuka, Ikeda, Hoshihira, Heinola, 
Liu and Ikeda-Ostuka, while hTST+TE and WZ+TE denote the calculated results based on the 


harmonic transition state theory and the Wert-Zener model with the temperature effect [30]. 


溶解 度 和 扩散 系数 是 多 中 氧 洁 留 研究 中 最 基本 的 物理 参数 ， 但 是 到 现在 为 止 ， 相 关 数 
据 十 分 稀少 且 存 在 较 大 差异 。 近 期 ， 我 们 采用 DFT 方法 研究 了 氧 的 扩散 和 溶解 性 质 ， 并 采 
用 准 简 谐 近似 方法 考察 了 温度 对 他 们 的 影响 [30]。 研 究 发 现 ， 温 度 无 法 改变 毛 在 多 中 的 最 
优 占 据 位 置 和 最 优 扩散 路 径 ， 但 是 会 显著 影响 溶解 能 和 扩散 激活 能 。 随 着 温度 的 增加 ， 和 氧 
x< Ze Sch EEE ELT KERTIS CRE ERM. EF DFT 计算 结果 ， 采 用 Sievert 定律 和 
= 过 渡 态 理论 计算 出 了 300 


2700K 内 氧 的 溶解 度 和 扩散 系数 。 我 们 的 计算 结果 在 高 温 范 围 内 (大 于 1500K) 与 实 
验 数 据 一 致 ， 在 低温 范围 内 与 实验 数据 有 较 大 差异 。 低 温 范围 内 的 差异 可 以 通过 缺陷 捕获 
效应 进行 解释 〈 图 4) 。 同 时 ， 我 们 发 现 与 未 进行 温度 修正 的 DFT 计算 结果 相 比 ， 温 度 修 
正 的 DFT 结果 与 实验 数据 符合 的 更 好 。 这 表明 在 采用 DFT 研究 钨 中 氧 滞留 性 质 时 ， 需 要 
对 计算 得 到 的 数据 进行 温度 修正 。 


图 5. (a)-(c): 间 际 氢 团 簇 结 构 示意 图 和 相应 结合 能 ;(d): 体系 吉 布 斯 自由 能 随 间 际 团 
大 小 的 变化 以 及 间隙 毛 团 簇 成 核 所 需 临 界 浓度 (e) 理论 预测 结果 与 实验 结果 对 比 〈 图 中 
数据 点 为 实验 值 ， 红 色 直 线 为 理论 预测 值 ) 。 


Figure 5. (a)-(c) configurations and binding energies of interstitial Hi。self-clusters. W atoms are 


in mineral blue and hydrogen atoms are in white. The red numbers represent spatial sequence of H 


occupation. (d) changes of the Gibbs free energy with the size of the interstitial H cluster. The red 
line highlights the critical concentration for H self-cluster nucleation. (e) comparison of the 
theoretical predictions with the experimental data (symbols denote experimental bubbling data and 


the red line represents the theoretical prediction). 
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氢 泡 的 成 核 机 制 却 不 清楚 。 近 期 ， 我 们 采用 理论 模拟 结合 热力 学 分 析 的 方法 研究 了 钨 基 材 
料 中 氧 自 团 复 行为 31]。 研 究 结果 表 明 间 隐 氧 团 复 行为 与 间隙 氨 团 簇 行为 类 似 ， 都 倾向 与 
形成 二 维 平面 板 状 结构 ， 且 间隙 氧 间 结 合 能 随 团 复 的 增 大 而 增 大 ， 最 终 收敛 到 一 个 常数 
(El Sac) 。 进 一 步 的 热力 学 分 析 发 现 由 于 间隙 氧 之 间 结 合作 用 较 弱 ， 通 常情 况 下 很 难 通 
过 自 捕获 形成 间 际 氧 团 簇 ， 但 是 在 氧 高 浓度 情况 下 ， 间 际 氧 可 以 聚集 形成 间隙 团 徐 (图 
5d) 。 这 一 发 现 为 低能 氧 缀 击 下 钨 基 材 料 中 氧 泡 形 成 提供 了 一 种 可 能 的 理论 解释 (图 Se)。 

近期 实验 结果 表明 辐 照会 在 钨 中 产生 的 大 量 缺 陷 ， 并 作为 氧 的 捕获 点 ， 致 使 氧 滞留 量 
显著 提高 ， 而 选择 合适 的 温度 和 合金 元 素 ( 如 Re、Ta、Ti 等 ) 能 抑制 辐 照 诱导 的 氢 滞 留 增 
加 。 因 此 ， 我 们 采用 DFT 研究 了 空位 -合金 元 素 - 氧 相 互 作用 参数 (结合 能 和 作用 半径 等 )、 
规律 以 及 合金 元 素 对 氢 滞 留 的 影响 [32,33]。 研 究 发 现 合 金 元 素 与 空位 、 自 间隙 都 有 较 强 的 
吸引 做 用 ， 其 中 与 自 间隙 作用 最 强 ， 可 以 阻碍 自 间隙 扩散 ， 合金 元 素 能 捕获 多 个 氧 原子 形 
成 间隙 氧 团 徐 ， 在 一 定 程度 上 降低 氧 原子 的 有 效 扩散 系数 ;合金 元 素 对 空位 - 氧 相 互 作用 影 
响 较 小 。 我 们 又 进一步 在 电子 层次 上 探讨 了 空位 -合金 元 素 、 自 间 隐 -合金 元 素 以 及 氧 -合金 
元 素 相互 作用 机 理 ， 并 对 作用 规律 开展 了 唯 象 讨论 。 发 现 ， 具有 较 大 电 负 性 的 溶质 原子 倾 
向 于 与 空位 结合 ， 有 具有 较 小 金属 半径 的 溶质 原子 倾向 于 与 自 间 隐 结合 ， 具有 较 大 原子 半径 
和 较 大 电 负 性 的 溶质 原子 与 氧 相互 作用 主要 受 化 学 相互 作用 控制 ， 其 结合 能 与 电 负 性 负 相 
关 ， 而 具有 较 小 原子 半径 和 较 大 电 负 性 的 溶质 原子 与 氧 相互 作用 主要 受 弹 性 相互 作用 控制 
其 结合 能 与 原子 半径 负 相 关 。 此 外 ， 我 们 采用 Simonovic 模型 和 经 典 Mac-Nabb and Forester 
方程 分 析 合 金 元 素 对 氧 分 布 和 扩散 的 影响 ， 并 结合 辐 照 点 缺陷 演化 行为 定性 探讨 了 合金 元 
素 对 点 缺陷 的 影响 。 发 现 合 金 元 素 一 方面 作为 固有 缺陷 ， 可 成 为 氨 的 滞留 点 ， 提 高 氧 滞留 
量 ， 另 一 方面 可 缩小 空位 和 自 间隙 扩散 速率 差 ， 促 进 辐 照 缺陷 复合 ， 降 低 辐 照 缺 陷 浓 度 ， 
进而 降低 氧 滞留 。 该 结果 定性 解释 了 钨 - 铂 、 钨 - 锂 合 金 中 氧 滞留 实验 结果 。 从 氧 滞留 、 合 
金 元 素 相 析出 和 巡 变 效应 的 角度 考虑 ， 我 们 认为 乌 是 综合 性 能 较 好 的 合金 元 素 。 


图 6. (a)-(c): 惰性 气体 间 际 团 
结构 示意 图 和 相应 结合 能 ，(d): 惰性 气体 富 集 层 阻 碍 氧 向 表 扩 散 和 释放 ， 增 大 氢 滞 留 示意 
图 [13,14]。 


Figure 6. (a)-(c) The stable configurations of inert gas atom clusters and their corresponding 


Ne 


binding energies in the simultaneous and sequential way; (d) schematics of the inert gas atoms 


enriched layer hindering H diffusion and release and increasing H retention [13, 14]. 
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杂质 He, Ne 和 Ar 倾向 于 沿 (110) 面 形成 稳定 的 二 维 平板 结构 (图 6a-b) 。 其 物理 机 理 是 间 
REI 


引起 局 域 晶 格 膨胀 ， 形 成 沿 (110) 面 扩展 的 低 电荷 密度 区 ， 有 利于 惰性 气体 溶解 和 聚集 
[34,35]。 惰 性 气体 团 复 随 聚集 原子 个 数 的 增 大 会 通过 激发 自 间隙 ， 突 变 成 更 加 稳定 的 空位 - 
惰性 气体 团 禾 。 由 于 惰性 气体 极 易 团 簇 化 并 形成 稳定 的 空位 -惰性 气体 团 徐 ， 因 此 在 聚变 堆 
中 惰性 气体 往往 会 在 钨 近 表 面 处 聚集 ， 形 成 惰性 气体 富 集 层 。 该 富 集 层 阻碍 氨 向 表面 扩散 
和 释放 ， 从 而 增加 了 和 氧 向 材料 内 部 的 扩散 ， 导 致 材料 内 部 氢 滞 留 量 增 加 《图 6c) 。 另 外 ， 
惰性 气体 占据 空位 后 会 降低 氢 在 空位 对 氢 的 捕获 能 力 ， 进 而 导致 辐 照 后 钨 表面 惰性 气体 富 
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co 相互 作用 参数 ， 我 们 构建 出 各 类 缺陷 的 能 量 参数 库 ， 并 输入 到 我 们 自主 开发 的 OKMC 
= 模拟 程序 中 。 如 图 7 所 示 ，OKMC 能 模拟 的 动力 学 过 程 包括 点 缺陷 的 复合 和 团 簇 化 ， 氧 在 
各 类 缺陷 附近 的 钉 扎 以 及 脱 钉 扎 ， 氧 、 点 缺陷 的 扩散 以 及 表面 托 脱 附 等 。 在 OKMC 模拟 中 ， 
© 程序 会 将 DET 计算 所 得 的 能 量 、 能 任 转 换 为 动力 学 速率 ， 从 而 以 该 速率 驱动 相应 的 动力 学 
一 过 程 。 在 本 工作 中 ， 我 们 用 SRM 等 蒙特 卡 罗 输 运 模拟 程序 ， 模 拟 了 高 能 所 离子 注入 纯 铭 
A 后 气 离子 的 深度 分 布 、 辐 照 缺 陷 的 深度 分 布 及 具体 缺陷 构 型 等 。 将 这 些 参数 输入 到 我 们 的 
> OKMC 模拟 程序 当中 ， 并 严格 按照 实验 条 件 设置 模拟 温度 、 辐 照 环境 、 时 间 等 变量 ， 从 而 


获得 能 与 实验 直接 对 比 的 模拟 结果 。 


<= Kl 7. OKMC 模型 及 模拟 中 涉及 的 动力 学 过 程 示意 图 。 


Figure 7. OKMC model and the dynamic processes involved in the simulations. 


首先 ， 我 们 按照 实验 文献 [36] 的 条 件 进 行 了 液 氨 温度 下 ， 低 能 气 (200 eV) 注入 后 的 
热 脱 附 模 拟 。 如 图 8 所 示 ， 模 拟 发 现 ， 随 着 温度 的 升 高 ， 间 隙 所 原子 的 迁移 在 100 K 左右 
被 激活 ， 对 应 第 一 个 脱 附 台阶 。 随 后 ， 在 200K. 400K 左右 ， 被 杂质 氧 原子 、 本 征 空位 钉 
扎 的 气 开始 挣脱 俘获 ， 分 别 对 应 第 三 和 三 个 台阶 。 模 拟 与 实验 结果 高 度 吻 合 。 随 后 我 们 模 
拟 了 更 为 贴近 聚变 堆 辐 照 环 境 的 实验 结果 ， 按 照 文中 [37] 中 的 实验 条 件 模 拟 气 原 子 从 表面 
脱 附 的 速率 。 我 们 模拟 得 到 的 结果 与 文献 [37] 的 实验 结果 吻合 很 好 ， 所 得 的 热 脱 附 峰 集中 
在 500K 左右 (如 图 9a 所 示 ) 。 此 外 ， 在 700K 附近 有 一 个 较 小 的 肩 峰 (用 蓝 色 竖 线 标 
H) 。 当 温度 升 至 800K 左右 ， 钨 中 和 气 原 子 已 经 完全 脱 附 出 去 。 同 时 ， 通 过 分 析 热 脱 附 谱 
及 缺陷 构 型 〈 如 图 9b 所 示 ) ， 我 们 发 现 500K 左右 的 低温 峰 是 由 气 原子 从 单 空位 中 脱 附 所 
引起 的 ， 而 TOOK 左右 的 高 温 峰 则 完全 由 空位 团 艇 所 贡献 。 
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图 8. 液 氨 温度 下 ， 少 量 200 eV 气 离 子 注入 后 ， 持 续 升温 时 钨 中 氧 同位 素 滞留 实验 值 与 
OKMC 模拟 值 对 比 。 


Figure 8. Comparison of experimental data of H isotope retention with the OKMC simulation 


results in W after 200 eV D* implantation at the liquid helium temperature. 


图 9. 室温 附近 下 ， 预 先 以 的 通 量 注 入 3 keV/200 eV 气 离子 至 剂量 后 ， 钨 中 和气 的 热 脱 附 谱 。 
(a) 模 拟 热 脱 附 谱 与 实验 结果 的 对 比 ; (b) 不 同 缺陷 类 型 对 脱 附 谱 的 贡献 。 


Figure 9. Thermodesorption spectra of D from W irradiated with 200 and 3 keV D*+ for low 
fluence of 10” D/m’ at a flux of 10” D/m’/s at room temperature. (a) comparison of simulation 
desorption sepectra with the experimental results; (b) decomposition of desorption spectra due to 


the contribution of different kinds of defects. 
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引起 钨 低温 脆性 的 主要 原因 之 一 是 杂质 元 素 偏 聚 在 晶 界 ， 造 成 晶 界 结合 强度 降低 ， 
致 晶 间 脆 断 [3, 38]。 而 微 合 金 化 和 晶 界 优化 设计 是 提高 材料 力学 性 能 的 重要 途径 ， 通 过 添 
加 一 些 有 益 的 合金 元 素来 强化 界面 、 降 低 界 面 能 ， 不 仅 可 有 效 提高 材料 的 强度 与 韧性 ， 而 
且 还 能 提高 晶 界 结构 的 热 稳定 性 。 然 而 ， 目 前 合金 元 素 强化 钨 唱 界 的 作用 规律 与 机 制 尚 未 
清楚 。 例 如 ， 哪 些 元 素 易 强化 / 脆 化 晶 界 ? 合金 元 素 强化 晶 界 能 力 与 唱 界 的 结构 与 类 型 是 否 
FERR? 因此 ， 人 们 需要 全 面 认识 合 金 元 素 强化 铝 唱 界 能 力 与 合金 元 素 自身 以 及 唱 界 结 
构 之 间 的 关系 ， 进 而 为 忽 基 材料 合金 元 素 的 优化 选择 以 及 蝇 界 结构 的 调控 提供 理论 指导 。 
基于 DEFT 方 法， 我 们 研究 了 19 种 过 渡 族 合金 元 素 (3d: Ti-Ni; 4d: Zr-Pd; 5d: Hf-Pt) 在 
不 同类 型 钨 对 称 倾 侧 晶 界 〈 倾 侧 轴 为 [100]，[110]，[111]) 中 的 偏 聚 和 强化 / 脆 化 效应 [39]。 
通过 能 量 学 方法 〈 偏 聚 能 和 强化 能 ) 和 动力 学 拉 伸 模拟 〔 拉 伸 强 度 ) 来 定量 描述 元 素 偏 聚 
引起 的 唱 界 结合 强度 的 变化 ， 获 得 强化 能 与 晶 界 结构 以 及 与 合金 元 素 自身 性 质 〈 人 金属 半径 / 
最 外 层 价 电子 数 ) 之 间 的 关系 。 研 究 结果 表明 合金 元 素 强化 / 脆 化 晶 界 能 力 与 晶 界 结构 密切 
相关 。 合 金 元 素 易 强 化 晶 界 能 较 大 的 晶 界 ， 而 易 脆 化 晶 界 能 较 小 的 晶 界 。 此 外 ， 合 金 元 素 
强化 唱 界 能 力 与 元 素 自身 的 金属 半径 成 正 相关 (如 图 10 所 示 ) 。 研 究 结 果 表 明 ， 尺 寸 效 应 
在 其 偏 聚 强化 界面 过 程 中 起 主导 作用 ， 即 金属 半径 比 忽 小 的 合金 元 素 ， 易 在 晶 界面 处 偏 聚 
并 能 强化 界面 。 该 研究 为 高 性 能 钨 基 材 料 合金 元 素 的 优化 选择 以 及 唱 界 的 优化 设计 提供 了 
理论 指导 : 在 晶 界 处 添加 Zr. Hf. Ta, Re 及 Ru 等 元 素 可 有 效 提 高 唱 界 的 结合 强度 ， 从 而 
改善 材料 的 力学 性 能 。 
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图 10. 不 同 唱 界 结构 中 合金 元 素 的 强化 能 与 元 素 金 属 半 径 及 唱 界 能 之 间 的 关系 图 [39]。 


Figure 10. Dependence of strengthening energies 4E sz of alloying solutes on the metallic radii of 


solutes as well as the grain boundary energy for different tungsten grain boundary structures [39]. 


一 般 来 说 ， 溶 质 偏 聚 导致 的 晶 界 结合 力 的 增强 /弱化 主要 有 以 下 两 种 机 制 : G) 尺寸 
效应 一 偏 聚 到 晶 界 附近 的 较 大 尺寸 的 原子 会 导致 晶 界 的 膨胀 ， 从 而 弱化 晶 界 结合 强度 ; 
Gi 电子 效应 一 通过 用 溶质 原子 取代 主体 原子 ， 原 有 的 主体 -主体 原子 间 结 合 键 将 被 更 强 
或 更 弱 的 主体 -溶质 原子 结合 键 所 取代 。 唱 界 脆 化 杂质 通常 将 相 邻 金属 原子 的 电荷 吸引 到 自 
身上 ， 削 弱 了 将 晶 界 保持 在 一 起 的 主体 -主体 结合 键 ; 而 唱 界 增强 杂质 则 不 会 从 金属 原子 上 
吸取 电荷 ， 反 而 它们 与 金属 原子 形成 相当 的 单 极 键 〈homopolar bond) ， 从 而 增强 唱 界 结 
合 强度 。 关 于 溶质 偏 聚 引起 的 晶 界 脆 化 /强化 的 机 制 ，Seah 很 早 提出 键 断裂 模型 (bond- 
breaking model) 来 预测 Fe 中 不 同 溶质 原子 的 强化 能 [40]。 后 来 ，Geng、Freeman 和 Olson 
提出 一 种 唯 象 模型 来 预测 不 同 金属 中 溶质 在 晶 界面 上 的 强化 能 [41]。 他 们 使 用 了 一 个 改进 
的 键 断裂 模型 ， 并 增加 了 弹性 失 配 项 ， 晶 通过 第 一 原理 计算 验证 了 模型 的 正确 性 [41]。 近 
年 来 ，Gibson 和 Schuh 系统 概述 了 不 同 金属 中 的 溶质 诱导 晶 界 结合 强度 的 变化 规律 ， 发 现 
溶质 诱导 唱 界 结合 强度 变化 的 主要 原因 可 以 通过 键 断裂 模型 来 解释 的 [42]。 最 近 ，Tran 等 
人 在 研究 溶质 原子 对 Mo 晶 界 结合 强度 变化 的 影响 规律 时 ， 提 出 两 参量 线性 模型 〈two- 
factor model) ， 即 强化 能 与 内 聚 能 变化 量 ( 全。)〉 和 人 金属 半径 相对 变化 量 ((Rx-R0)/Ro) 之 
间 成 线性 关系 [43]。 他 们 发 现 两 参量 线性 模型 可 以 很 好 描述 不 同 溶质 原子 在 Mo 唱 界 处 的 强 
化 能 。 通 过 综合 分 析 W, Mo, Ni, Fe 等 唱 界 处 合金 元 素 的 强化 能 的 变化 规律 ， 我 们 发 现 
线性 模型 可 以 拓展 应 用 于 描述 不 同 金属 晶 界 处 合金 元 素 的 强化 能 (如 图 11 所 示 ) 。 
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CS 图 11. 不 同 合金 原子 在 23(111) W 唱 界 处 偏 聚 的 强化 能 计算 值 与 两 参量 模型 拟 合 值 对 比 。 


Figure 11. Comparison of the calculated values of strengthening energies of different solute atoms 


at the £3(111) W grain boundary with the predictions by two-parameter model. 
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聚变 环境 下 ， 高 能 中 子 〈14MeV) 辐 照 往往 导致 W 中 产生 Cn, y) 和 Cn, 2n) We 
反应 而 产生 Ta, Re 和 Os 等 同位 素 原子 。 最 近 理 论 研 究 发 现 ， 在 未 来 聚变 环境 下 服役 五 年 ， 
£5 HW AS 7G Ta, Re 和 Os 浓度 分 别 达 0.84at% 、3.8at% 和 1.4at%[44]。 中 子 辐 照 又 导致 妨 
变 元 素 Ta、Re 和 Os 在 钨 中 发 生 聚 集 而 形成 偏 聚 物 。Nemoto 等 人 和 Tanno 等 人 [45-47] 研 究 
发 现 ，W-Re 二 元 合金 在 高 能 中 子 辐 照 下 材料 中 出 现 盘 状 和 条 状 两 种 Re 偏 聚 物 ， 而 偏 聚 物 
又 导致 钨 材料 的 硬化 和 脐 化 。 然 而 ， 氢 / 氛 气 泡 及 偏 聚 物 形 成 的 微观 机 理 尚 不 清楚 。 此 外 ， 
巡 变 元 素 和 氧 / 氨 共 同 存在 于 面向 等 离子 体 材料 中 ， 二 者 之 间 相 互 作用 机 制 也 不 清楚 。 针 对 
这 些 问题 ， 我 们 使 用 DFT、 集 团 展开 及 动力 学 蒙特 卡 罗 等 方法 研究 了 钨 中 氧 / 氨 气泡 、 娘 变 


元 素 偏 聚 物 形成 及 二 者 之 间 的 相互 微观 机 理 [48-50]。 
基于 DFT 方法 ， 我 们 研究 了 钨 中 氧 / 氨 气 光 成 核 和 生长 机 制 [48]。 研 究 发 现 ， 氧 原子 很 
难 在 W 的 间 队 位置 聚 集 ， 即 当 钨 中 注入 氧 浓度 较 低 时 〈 即 小 于 临界 浓度 ) ， 此 时 精 效 应 起 
主导 作用 (如 图 5c-d 所 示 ) 。 由 于 温度 的 影响 ， 由 氢 原 子 组 成 的 小 团 筷 很 容易 瓦解 。 当 所 
浓度 较 大 于 临界 浓度 时 ， 业 效 应 对 小 团 艇 仍 起 主要 作用 ;但 对 较 大 团 和 伐 ， 此 时 答 起 主导 作 
用 ， 导 致 吉 布 斯 自由 能 随 氮 原 子 个 数 的 增 大 呈现 先 增加 后 下 降 ， 使 得 体系 中 较 大 毛 团 簇 很 
容易 聚集 ， 这 些 氧 团 马 可 作为 纯 钨 中 氧 泡 的 成 核 点 。 对 于 氨 原 子 ， 它 们 倾向 于 在 铭 的 
(110) 层 间 聚 集 形成 密 堆 积 的 单 层 结构 。 单 层 中 氨 原 子 之 间 最 小 距离 约 1.6A。 氨 单 层 的 
结合 能 随 着 氨 原 子 数 指数 增加 ， 当 氨 原 子 数 为 9 时， 结合 能 高 达 16 eV。 氨 原子 之 间 几乎 无 
E 子 相互 作用 。 氢 单 层 可 能 为 实验 上 观察 到 的 扁平 状 氨 气 泡 的 初始 结构 。 由 于 局 域 应 力 的 

存在 ， 当 氨 单 层 中 的 氨 原 子 数 大 于 5 时， 氨 单 层 的 附近 可 能 发 射 自 间 阶 原子 而 形成 弗兰克 
一 尔 缺 陷 对 。 氨 原子 则 可 以 自发 地 扩散 到 该 缺陷 对 中 的 位 形成 空位 - 氨 复 合体 。 此 外 ， 研 究 发 
> 现 所 和 和 氮 均 倾向 于 在 空位 中 占据 。 基 于 质量 反应 定律 ， 我 们 计算 了 空位 -所 复合 体 的 浓度 随 
温度 的 变化 关系 。 该 结果 表明 在 ITER 的 服役 温度 (1000K) 附近 ， 空 位 - 氧 复合 体 稳定 存 
在 ， 尤 其 是 空位 捕获 不 多 于 6 个 氧 的 复合 体 大 量 存在 ， 而 空位 - 氨 复 合体 则 更 加 稳定 。 我 们 
研究 还 表明 ， 氢 和 氨 在 空位 中 的 聚集 增 大 了 空位 与 空位 之 间 的 结合 能 。 当 空位 中 聚集 的 所 
原子 数 大 于 9 时 ， 该 空位 近邻 可 能 发 射 一 个 自 间隙 原子 。 因 此 ， 氢 和 氨 气 泡 可 以 通过 以 下 
两 种 方式 成 核 和 生长 : (1) 空位 捕获 所 原子 ， 氨 原子 浓度 达到 阔 值 时 ， 该 空位 周围 发 射 一 
个 自 间隙 原子 或 ， 长 大 的 空位 再 捕获 氨 原 子 再 发 射 自 间隙 原子 直到 氮气 泡 的 形成 。 (2) A 
和 和 氨 原子 在 空位 中 占据 ， 促 进 了 空位 的 聚集 从 而 形成 空位 团 秘 ， 而 氢 和 氨 在 该 空位 团 簇 中 
的 聚集 导致 气 和 所 气泡 的 形成 。 


PT 图 12. W 中 Ta, Re 和 Os 在 无 和 有 空位 情况 下 偏 聚 物 结合 能 随 原 子 数 变化 。 


Figure 12. Binding energies of (a) Tan, Ren, and Os, clusters without vacancy and (b) VacTan, 


VacRe,, and VacOs, with a vacancy as a function of the number of solute atoms in W. 


此 外 ， 我 们 使 用 DFT 方法 和 集团 展开 方法 系统 研究 了 W-Ta, W-Re 及 W-Os 等 合金 中 
妨 变 元 素 偏 聚 物 形成 的 微观 机 理 以 及 Os 和 Ta 等 对 Re 偏 聚 物 形成 的 影响 [50]。 研 究 发 现 ， 
无 辐 照 空位 时 仅 Os TE HE, TT Ta 和 Re 均 无 法 偏 聚 。 辐 照 空位 的 存在 促进 了 娘 变 元 
WEE, Re 和 0s 均 可 在 有 辐 照 空位 的 钨 中 偏 聚 ， 而 Ta 则 无 法 偏 聚 。 单 空位 最 多 可 以 结 
合 8 个 Re 原子 和 至 少 14 个 Os 原子 (如 图 12 所 示 ) 。 空 位 -Os 的 结合 强度 要 远大 于 空位 - 
Re 的 结合 强度 。 这 些 研究 结果 解释 了 实验 发 现 的 在 中 子 辐 照 条 件 下 ，Re TE HMR, 
而 Ta 不 在 钨 中 偏 聚 的 现象 。 我 们 的 研究 还 发 现 ， 巡 变 元 素 Os 可 促进 Ta 和 Re 的 偏 聚 ， 而 
Ta 的 存在 可 抑制 Re 和 Os 的 偏 聚 ， 这 一 研究 结果 也 与 实验 观察 到 的 Ta 可 抑制 Re 的 偏 聚 等 


结果 定性 一 致 。Os A ED ATR EB ti ET res MOSS, EE 500-700K 之 间 时 ， 达 到 
最 大 值 。 


图 13. Ta, W, Rei Os 中 He 团 徐 的 (a) 总 结合 能 和 (b) 平 均 结合 能 随 He 聚集 个 数 之 间 的 变 
化 关系 。(a) 中 空心 符号 是 考虑 了 零点 能 量 修正 后 的 总 结合 能 。 


Figure 13. (a) Total binding and (b) averaged energies of He, clusters in Ta, W, Re and Os as a 
function of the number of solute atoms. The solid and hollow symbols in (a) are total binding 
energies without and with ZPE corrections, respectively. 


在 此 基础 上 ， 我 们 还 研究 了 钨 中 偏 聚 物 对 氮气 泡 形成 的 影响 。 发 现 ， 所 在 Ta 偏 聚 物 中 
占据 八 面体 位 置 最 稳定 ， 且 其 最 优 扩 散 路 径 是 从 一 个 四 面体 到 最 近邻 四 面体 位 置 。 与 在 铅 
中 的 行为 类 似 ， 氨 在 Ta 的 (110〉 层 间 间 际 位 置 察 集 形成 密 堆 积 单 层 结构 。 氨 则 倾向 于 占 
据 Re 和 Os 偏 聚 物 的 基底 八 面体 位 置 ， 其 最 佳 移动 路 径 是 从 基底 八 面 体位 置 沿 <0001> 方 向 
经 过 八 面体 位 置 扩散 到 最 近邻 基底 八 面 体位 置 。 氨 在 Re 和 Os 的 间隙 位 置 聚集 成 核 形成 无 
对 称 构 型 的 氨 团 禾 。 氨 在 Ta、Re 和 Os 等 侦 聚 物 中 聚集 降低 所 附近 格 点 原子 的 稳定 性 而 发 
射 格 点 原子 从 而 导致 氨 团 簇 长 大 成 泡 。 经 过 比较 发 现 ， 铅 中 氨 的 间 际 团 

结合 强度 最 大 ， 而 Os. Re 及 Ta 等 偏 聚 物 中 氨 间 隙 团 

结合 强度 相对 较 弱 ， 且 Ta Pel HI 

结合 强度 最 弱 〈 如 图 13a) 。 辐 照 空位 的 存在 促进 了 氢 在 合金 元 素 偏 聚 物 中 的 聚集 ， 
但 是 空位 - 氛 团 簇 的 结合 强度 顺序 依然 是 W>Os>Re>Ta〈 如 13b) 。 与 空位 类 似 ， 自 间隙 原 
子 也 可 成 为 氨 原 子 的 捕获 点 。 而 当 Ta、Re 和 Os 与 自 间隙 原子 结合 形成 混合 自 间 隙 原子 时 ， 
氨 与 混合 自 间隙 的 结合 强度 要 低 于 氮 与 自 间 隙 的 结合 强度 ， 且 所 与 Ta 混合 自 间隙 的 结合 强 
度 最 小 。 因 此 ， 钨 中 Ta, Re 和 Os 等 偏 聚 物 的 存在 一 定 程度 上 降低 了 氢 的 结合 强度 ， 从 而 
al = ME TUE FRR BOA ATR KE 


六 、 结 束 语 

材料 的 辐 照 损伤 是 一 个 多 尺度 模拟 问题 。 它 们 涉及 从 原子 尺度 到 连续 介质 尺度 上 、 从 
纳 秒 到 秒 时 间 尺 度 上 不 同 物理 机 制 的 耦合 和 关联 。 咪 需 运 用 多 尺度 模拟 来 获取 材料 在 辐 照 
环境 下 微观 、 介 观 和 宏观 等 不 同 层次 的 结构 演化 信息 。 近 十 多 年 来 ， 随 着 并 行 计算 科学 与 
技术 的 迅速 发 展 ， 多 尺度 模拟 方法 已 成 为 核 聚 变 堆 中 等 离子 体 与 壁 材 料 相互 作用 研究 领域 
的 一 个 可 靠 的 理论 手段 ， 并 起 到 了 不 可 或 缺 的 重要 作用 。 国 内 聚变 堆 材料 计算 模拟 工作 ， 
近年 来 特别 是 ITER 专项 设置 以 来 ， 无 论 从 计算 模拟 研究 还 是 新 的 计算 方法 开发 应 用 方面 
都 取得 了 长 足 的 进步 和 发 展 。 本 文 主要 介绍 了 近期 我 们 课题 组 在 多 尺度 模拟 钨 基 材 料 辐 昭 
损伤 与 氧 氨 效应 机 理 方面 的 一 些 进展 ， 主 要 包括 纳米 唱 铝 / 铁 的 抗 辐 照 损 伤 与 自修 复 机 理 、 
钨 中 氢 清 留 机 制 、 合 金 元 素 强化 界面 机 制 、 妨 变 元 素 仿 聚 及 其 对 氢 泡 形成 机 制 的 影响 等 内 


容 。 通 过 近 八 年 的 研究 ， 我 们 自主 开发 了 一 套 多 尺度 模拟 程序 ， 搭 建 了 模拟 核 材 料 辐 照 


2i 


伤 的 软件 平台 〈 温 度 加 速 分 子 动力 学 、 静 态 计算 、 基 于 第 一 性 原理 的 集团 以 及 实体 动力 学 
蒙特 卡 洛 方法 等 ) ， 基 于 该 模拟 平台 ， 我 们 揭示 了 纳米 晶 材料 中 晶 界 附近 原子 通过 “ 巡 
游 ” 来 修复 辐 照 损伤 的 自修 复 机 制 ， 为 研发 高 性 能 纳米 晶 抗 结 照 材料 提供 了 新 思路 ， 建 立 
了 钨 中 氧 / 氨 - 空 位 -合金 / 妨 变 元 素 相互 作用 数据 库 ， 利 用 OKMC PEPFAR AR TARRA 


热 脱 附 谱 的 影响 机 理 ; 构建 了 钨 品 界 合金 元 素 强化 界面 与 唱 界 类 型 以 及 合金 元 素 自 身 性 质 
之 间 的 定量 关系 ， 为 高 性 能 龟 基 材料 晶 界 的 优化 设计 与 合金 元 素 的 优化 选择 提供 了 理论 指 
导 。 然 而 ， 必 须 指出 ， 利 用 计算 模拟 有 效 获得 14MeV 中 子 辐 照 下 材料 结构 演化 信息 仍然 是 
一 项 巨大 的 挑战 。 此 外 ， 未 来 我 们 需要 进一步 实现 不 同 尺度 之 间 模 拟 方法 的 有 效 衔 接 以 及 
关注 多 场 耦合 作用 下 的 材料 服役 性 能 评价 及 其 演化 规律 研究 ， 从 而 理解 并 揭示 材料 的 辐 昭 
损伤 机 理 及 性 能 失效 机 制 ， 为 中 国 聚 变 工程 实验 堆 第 一 壁 选材 提供 科学 依据 。 
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Recent Progress on Mechanisms of Irradiation Damage and Hydrogen/Helium 
Behaviors of Plasma-Facing Tungsten in Nuclear Fusion Reactors by Modeling 
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Abstract: Tungsten-based materials have been considered as one of the most promising 
candidates for plasma facing materials in nuclear fusion reactors because of the high melting 
temperature, high thermal conductivity, good resistance to neutron irradiation and high sputtering 
resistance. In a fusion environment, the bombardment of high energy (14 MeV) neutron and high- 
flux and low-energy hydrogen (H) and helium (He) ions leads to severe radiation damage of 
tungsten materials. Investigation on the mechanisms of irradiation damage and H/He behaviors in 
materials is significantly important to reveal the changes of microstructure and properties due to 
irradiation and also important to develop advanced irradiation-resistant materials. In recent years, 
with the development of computer simulation technology, multi-scale simulation techniques have 
been widely used to study the mechanism of irradiation damage and H/He behaviors of key 
materials in nuclear fusion reactors. In this work, in view of our recent years of experience, we 
mainly introduced the progresses on mechanisms of radiation-tolerance of nanocrystalline iron 
(Fe) and tungsten (W), H retention in W, strengthening effect of alloying elements on W grain 
boundaries, segregation of transmutation elements and their effect on He bubble formation in W. 
The structural evolution information of materials at different time/space levels under an irradiation 
environment has been obtained by multi-scale simulation, and so the mechanisms of irradiation 
damage and H/He effects were revealed, which provides theoretical guidance for the optimal 


design and performance control of plasma-facing material in nuclear fusion reactors. 
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helium behaviors; Multi-scale modeling and simulation. 


